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SUMMARY 

This report summarizes a study effort directed toward developing 
handbook- type methods for predicting, without undue conservatism, aero- 
dynamic heating rates to cones and cylinders at angle of attack. Applica- 
bility of these isolated body methods to composite body geometries consisting 
of cones and cylinders is also discussed. Approximate methods for compres- 
sion and expansion corner regions are also given. 

In general, the test data trends were predicted well by several theories. 
Laminar convective heating rates to cylinders at zero-degree angle of attack 
were well predicted by flat plate theory. For cylinders at angle of attack, 
yawed infinite cylinder theory applied. The laminar stagnation line heating 
rates were well predicted by the method of Kemp, Rose and Detra at the 
angles of attack less than 30 degrees. The turbulent stagnation line heating 
rates were well predicted by the method of Beckwith and Gallagher. For 
cones at angle of attack, the laminar and turbulent heating rates were well 
predicted by streamline divergence theory. Heating rate distributions off 
the stagnation line for cones and cylinders are also presented based on avail- 
able experimental data. Then, using the heat transfer theories that corre- 
lated best with the test data, extrapolations to typical boost trajectory flight 
conditions were made and the results are presented in the form of graphs 
that can be used for rapid estimates of aerodynamic heating rates. To aid 
in the calculation of aerodynamic heating rates, the local flowfield prpperties 
for cones and cylinders are also presented. Finally two example calculations 
are made to illustrate the use of the graphs presented. 
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NOMENCLATURE 


Symbols 


CP 

C 

P . 

D 

E 



pressure coefficient 
specific beat 

diameter 

crossflow momentum thickness ratio 
defined by Eq. (2.7) 


E, 


g 

H 

h 

h 

c 

k 

M 

M . 
wt 

N 

Nu 


P 

Pr 


q 

q 

R 

r,x 

Re 

St 

T 



defined by Eq, (2,18) 

gravitational constant 
total enthalpy 
static enthalpy 
heat transfer coefficient 

thermal conductivity 
Mach number 
molecular weight 

defined by Eq. (2.5) 

Nusselt number 
pressure 
Pranti number 

velocity used in Section 2.3.2 
convection heating rate 
radius of cylinder 

position coordinates used in Section 2.3,2 

Reynolds number 

Stanton number 

temperature (^R) 

recovery temperature 

normal component of freestream velocity 
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NOMENCLATURE (Continued) 

Symbols 

V velocity 

tangential component of freestream velocity 

X-. transformed characteristic length (laminar flow) 

transformed characteristic length (turbulent flow) 

X running length; measured along surface 

Z compressibility factor 


Greek 

a 

r 

6, p 

A 

e 

e 

p 

<t> 

e 

Subscripts 

G 

d 

D 

e 

Lr 

m, c 

o 

R 

r 

s 


angle of attack 

ratio of specific heats 

defined in Section 2.3.2 

sin ^ ( -^ ) , Mach angle 

attached shock angle 
flow angle 
density 

angle measured perpendicular to the body axis 
((t> is zero at the stagnation line) 

cone half- angle 

crossflow pressure gradient effects 
downstream of corner 

based on diameter 
boundary layer edge values 
local conditions 

evaluation at local pressure and mean boundary 
layer enthalpy 

evaluated at zero Mach number 
recovery value 

Hanks Rho Mu Reference condition (simplified method) 
stagnation conditions 

vi 
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Subscripts 


Superscript 


NOMENCLiATURE (Continued) 


stagnation line 
upstream of corner 
wall conditions 
based on x 

free stream condition ahead of shock 
right characteristic 
left characteristic 
property upstream of shock 
property downstream of shock 


Eckert Reference Enthalpy condition 
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Section 1 
INTRODUCTION 

Plans are being considered to launch several Saturn vehicles into polar 
orbit from the East Coast. This calls for launch maneuvers which require 
the vehicle to fly at angles of attack of up to approximately 50 degrees. 
Because Saturn vehicles were not originally analyzed for the thermal environ- 
ment created by these maneuvers, a reevaluation of the ascent phase heating 
must be made. There are at present, however, no completely effective means 
to predict induced aerodynamic heating at these large angles of attack. As a 
result, extremely conservative assumptions must be made which will result 
in excessive weight and cost penalties. It is therefore expedient that tech- 
niques be developed to predict these heating rates as accurately as possible 
at large angles of attack. 

In addition to flying the Saturn vehicle, NASA is currently involved 
with the development of a reusable Space Shuttle vehicle. The results of the 
efforts presented herein are also directly applicable to determining heating to 
various areas of the Shuttle vehicle. These areas include: (1) swept wing 
leading edges; (2) swept fin leading edges; (3) centerline of round bottom delta 
wings at angle of attack; and (4) bottom centerline of cylindrical fuselages. 

The following sections present a technical discussion of the work done 
under this contract. First, a discussion of the acquisition, reduction, and 
correlation of the test data is given. Next, results of the reductions and cor- 
relations are given. This is followed wdth the design curves for application 
of these results to flight heating calculations. 
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Section Z 

TECHNICAL DISCUSSION 

Discussed in this section are the sources of test data and methods used 
in the data correlations, the theoretical methods used for comparision with 
test data, and the methods used to calculate the local flowfield properties. 

It is evident that many sources of data exist in the literature. However, 
when conducting data correlations one has to be very selective in the choice 
of data used in order to isolate the various effects on heat transfer such as 
Reynolds number, Mach number, nose bluntness, angle of attack, wall tem- 
perature, compression and expansion corners, and other pressure gradient 
effects such as might exist in some test facilities. One also has to look at 
the accuracy of the instrumentation used to obtain the test data. 

To provide meaningful insight into the effects of some of the above vari- 
ables on heat transfer, several convective heat transfer theories were inves- 
tigated, Some of the theories were programmed and their results were com- 
pared with the data. The theories that correlated best were then used to 
provide extrapolations of the existing data to provide the "best" estimate of 
the aerodynamic heating environment expected during a typical boost trajec- 
tory. The theories were applied only to the stagnation line. The heating 
rate distributions off the stagnation line were obtained from actual test data 
and as such are limited in their applicability to conditions other than the test 
conditions. 

The local flowfield properties were needed to perform the data cor- 
relations. The method used in the literature for presenting the test data in 
many cases varied from source to source even for the same type bodies. It 
was decided in the course of this work to use two standard methods, one for 
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cones and another for cylinders. The flowfield methods used in this work 
are documented in Section 2.3 and flowfield results for a variety of flight 
conditions are presented in Section 4, 

2.1 ACQUISITION OF TEST DATA 

Data were obtained from numerous sources (Refs. 1 through 11). Cylin- 
der test data correlations were made using both local and freestream flow 
properties for yawed cylinders. The method used to obtain the local flow 
properties for cylinders is discussed in Section 2.3. Since there was no ap- 
parent gain in accuracy and elimination of data scatter to be gained using 
local flow properties all the cylinder data> correlations are presented, in 
Section 3.1 using freestream properties. 

Data correlations for cones were also reduced using both local and 
freestream properties. Since much less data scatter was evident using the 
local flow properties, the results are presented in this form. For cones at 
angle of attack the equivalent cone concept was used to obtain local flow 
properties along the most windward (stagnation) streamline. The numerical 
solution used for obtaining the conical flow properties is presented in Section 
2.3. 


Numerous other reports contained test data that were reviewed but 
were rejected for various reasons. Some of the reasons for rejecting them 
were: (1) insufficient data on test conditions; (2) uncertainty in methods used 

to obtain local flow properties in reducing data presentation; (3) abscence of 
other pertinent parameters needed for data correlations; (4) uncertainty in 
air properties used, i.e., conductivity and viscosity; (5) data were not appli- 
cable (not isolated body data); and (6) data were not clearly laminar or turbu- 
lent, but apparently transitional. 
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2,2 THEORETICAL, CONVECTIVE HEAT TRANSFER METHODS 

The boundary layer on axisymmetric bodies at angle of attack is three- 
dimensional. Solution of the general three-dimensional boundary layer equa- 
tions was not within the scope of this contract and, therefore, approximate 
techniques were used to provide first-order corrections for the effects of 
crossflow. 

2.2.1 Cones 

Several methods were tried and the method ultimately selected to cal- 
culate the convective heating rates along the most windv/ard streamline on 
cones was based on Streamline Divergence theory as first proposed by Vaglio- 
Laurin in Refs. 12 and 13. Since Vaglio-L»aurin* s initial work, many engin- 
neers have used the Streamline Divergence method with numerous modifica- 
tions. The method used here is still another modification. The modification 
was necessary only for the laminar boundary layer case in order to obtain 
good correlations at angle of attack. 

The basic idea of the Streamline Divergence method is relatively 
straightforward and is a natural continuation of the axisymmetric flat plate 
analysis. As is known, axisymmetric bodies at zero-degree angle of attack 
have been analyzed successfully by two-dimensional methods with the intro- 
duction of various characteristic length transformations. In these trans- 
formations, the local normal radius of curvature of the body is the major 
contributing factor to the *'stretched” characteristic length used in the heating 
equations. This arises from the momentum change in the boundary layer due 
to the spreading effect as it moves around the axisymmetric body. When the 
non- axisymmetric body is considered, or vehicle at an angle of attack in this 
case, the momentum change is no longer about the normal curvature of the 
axisymmetric surface, but is along the curvature of the local fluid stream- 
line. The corresponding length transformation was calculated in this work 
for both laminar and turbulent flow along the stagnation (most windward) 
streamlines. 
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The Streamline Divergence method used to evaluate this transform is 
explained in the following sections. 


Laminar Cones 


The basic task of the Streamline Divergence method is to provide a 
transformed characteristic length along the streamlines which will modify 
the familiar zero pressure gradient (flat plate) heating rate equations into 
relationships that will handle flows with both circumferential and longitud- 
inal pressure gradients. To do this we start with the following equation 
for zero pressure gradient flow fields: 


q 


-2/3 

0.332 (Hj. - h^) Pr* 



, * *.1/2 
(PM) 


( 2 . 1 ) 


This is the Blasius incompressible solution modified for compressible flow 
by the introduction of Eckert reference condition fluid properties, denoted 
by the asterisk {*). 

In addition, for comparison purposes, the following flat plate heating 
rate equation was also used, 


q = 0.332 (H - h ) 

r w' r 


ft)" 


(Pr 


1/2 


( 2 . 2 ) 


This again is the Blasius incompressible solution modified for compressible 
flow by the introduction of reference conditions based on Hanks Rho- Mu 
method (Ref. 14). 
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For an unyawed sharp cone in laminar flow, the solution for the char- 
acteristic length Xt is 

i-; 


X := 
L, 



(2.3) 


For the stagnation line of a yawed sharp cone (negligible streamwise 
pressure gradients), the solution for the transformed characteristic length 
is 


X = — (2.4) 

l + 2(l + NEj_^) 

where the parameters N and are calculated as follows (Ref. 14). The 
solution for the parameter N was obtained from the solution of the cross- 
flow momentum equation for inviscid flow along with certain simplifying 
assumptions which appear valid along the stagnation line. The result is. 




whe re 





(2.5) 


( 2 . 6 ) 



The quantity E_ is defined as 


= 1 + 0.718 1)(2 (2.7) 
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where 


exp K = 0 when N < 0.05 and 0.99 < N _< 1.01 

exp K = “0. 194 exp j - Z/3 N (N - 1) j when 0.05 < N <0,99 

exp K = 0.194 exp |- Z/3 (N - 1) j when N > 1.01 


and 


F 


S, c 


0.355 

0.402 r 


( 2 . 8 ) 


(ZT) 


m, c 


(ZT) 


e f SXj 


(2.9) 


The quantity (ZT)^ is evaluated at the local pressure and mean boundary 

f c 

layer enthalpy. Equations (2.7) through (2.9) were developed in Ref. 14 on 
the basis of providing the best fit to exact similiar solutions. 


# Turbulent Yawed Cones 

For turbulent flow over cones, three heating rate methods were 
selected for comparison with test data. They are: (1) the Eckert Reference 
Enthalpy method using the Blasius skin friction law, i.e., 


. ^ 0.0296 

^ . 0 . 666 ? 




Pr 





(H.-V 


( 2 . 10 ) 


(2) the Eckert Reference Enthalpy method using the Shultz - Grunow skin friction 
law, i.e., 



0.185 p* 

*0.6667 f * , 2.584 

Pr [Alog(Re )| 


( 2 . 11 ) 
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where 



* 

M 


( 2 . 12 ) 


and finally (3) the Hanks Rho-Mu method, also using the Shultz - Grunow skin 
friction law, i.e., 


q = 


0.185 Pr Ye 


_ 0.6667 FT! TZ , Qnnn\l 2.584 
Pr^^ IAlog(Re^ + 3000) 




^ ) 

w' 


(2.13) 


where, 




and 



Pr Mr^eX. 



(2.14) 



(2.15) 


For the stagnation line of an unyav/ed sharp cone in turbulent flow, the 
solution for the transformed characteristic length Xrp is, . 

= (|-) X (2.16) 
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For the stagnation line of a yawed sharp cone in turbulent flow, the 
solution for the transformed characteristic length is 


X s H 2 (2.17) 

^ l+-|U+E^N) 

where N is evaluated using Eq, (2.5), The parameter E^ is evaluated in Ref. 14 
as follows, 


E = 1 + 0.55 (l[ 


1 + F, 


S, c, o 




exp K 


E. 


T4 


1 + 0.718 (xTITfI - 1) (2S )6xp K 

^ S, C, O ' ' i- 


C, O' 


(2.18) 


where the quantities F^ ^ ^ and ^ are evaluated using Eqs. (2.8) and (2.9), 

respectively, except that the Mach number is now zero, therefore the mean 
boundary layer enthalpy is now 


h 


m, o 


0.5(h +h ) 
' e w' 


2.2.2 Cylinders 


(2.19) 


The methods used to calculate the three-dimensional boundary layer 
heat transfer on the stagnation line of yawed infinite cylinders are based on the 
Crossflow Approximation method. This method employs existing similiar 
solutions to the boundary layer equation. 
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• Laminar Yawed Cylinders 

Two stagnation line laminar heat transfer methods were programmed 
for comparison with the test data. One method used the correlations of Kemp, 
Rose and Detra (Ref. 15), and the other used Beckwith's yawed cylinder method 
(Ref, 16). The stagnation line convective heating rates fo,r equilibrium air 
was analyzed in Ref. 16 . The resulting heat transfer correlation obtained 
was 



0,565 


Pr 




w 




(H^-h ) 
' R w' 


( 2 , 20 ) 


If the crossflow velocity component is supersonic, then the velocity gradient 
is computed from Newtonian theory, i.e,, 


^ - JL 

dR ” R 



( 2 . 21 ) 


If the crossflow velocity component is subsonic, then the velocity gradient is 
computed by the following equation from Ref. 16: 


A^r sin a 

dv _ (30 

dR ^ 2R 


5.1 


cos(M sina) - 

(30 ' 



( 2 . 22 ) 


The other laminar heat transfer solution for the stagnation line of a yawed 
infinite cylinder is that of Beckwith (Ref. 16), i.e., 
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q 


0.354 



dV 

f sin a 

dR 



(’■r- ‘w ' 


(2 


Equation (2.23) is based on exact boundary layer solutions for constant spe- 
cific heat and a constant Prandtl number equal to 0.7. The velocity gradient 
is calculated using either Eq. (2.21) or (2.22). 


• Turbulent Yawed Cylinders 

The turbulent yawed cylinder stagnation line heat transfer method used 
was that presented by Beckwith and Gallagher in Ref. 8, i.e.. 


4 = 0.0288 


( V cosa\°*^ /j,r\0.2 « Q k 

-22 1 tp n -2° (T - 

U. / \dR) n 


T ) 

w' 


(2 


2.3 LOCAL FLOWFIELD CALCULATION METHODS 

This section presents methods used to calculate the local flowfield 
properties on cones and cylinders. Inviscid analysis methods were used, 
and no boundary layer interactions were considered. 

2.3.1 Yawed Cylinder Flow Fields 

The local flowfield at the stagnation line of a yawed infinite cylinder 
is computed as the real gas stagnation conditions behind a normal shock 
wave. The real gas (equilibrium air) properties are computed using the 
method of Hansen (Ref. 17). The normal shock calculation used is simillar 
to that presented in Ref. 18. Basically the solution involves the solution 
of the continuity, momentum, and energy equations plus an equation of state 
which includes compressibility effects. The following sketch shows the flow 
field geometry. 
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The pertinent equations are 


2 

IP 2 2 I Conservation 

^ ^ ^00 P ' \~/ Momenttim 



C ons e rvation 
of Energy 


Pi Ml - P2 M2 


C ontinuity 



Pz '^z ^2 


Equation of 
State 


2-11 


LOCKHEED vHUNTSVILLE RESEARCH & ENGINEERING CENTER 




L.MSC/HREC D162315 


An iteration procedure is used to obtain the solution to the preceding set of 
equations. After the normal shock solution is obtained, the flow is isentrop- 
ically decelerated such that u^ equals zero with a corresponding increase in 
and 

2.3.2 Cone Flow Fields 

The method used to solve the conical flow field is the method of Ref. 19. 
For cones at angle of attack the equivalent cone concept was used to obtain 
conical flow solutions along the most windward streamline. The method is 
well suited to treat ideal, frozen, or equilibrium reacting gas mixtures. The 
resultant equations require a double -iterative solution rather than an iterative - 
integral solution such as the Taylor- Mac coll method* 

Consider a right circular cone immersed in a supersonic flow (Fig. 2-1). 
The flow downstream of the shock wave is irrotational for conical flow. If 
the Mach number downstream of the shock wave is greater than unity, the 
characteristic equations are real. They are: 


dr 

dx 


tan(^ + A) 


(2.30) 


d(e + cotA^ + (^ = 0 

sin (I + A) ' 


(2.31) 


Along the right (-) characteristic the finite difference analog of Eq. (2.32) 
is written; 


t, - =0 

2 ^ (I A) Vv- 


sin 


(2.32) 
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while along the left (+) characteristic 



0 - cot A 


<^12 ■ 


sin A sin ^ 
sin(g - A) 



(2.33) 


where the barred values are average between the shock and the cone surface. 
With the aid of the physical characteristic equation (2.30) and geometry the 
axisymmetric term is found; 



(2.34) 


whe re 6 = 


cos (e + 0) sincr 
cos (a+ p) sine 


and where 


p -f = tan 


I cot(g + A) 


The final relations necessary for the solution are provided by oblique 
shock theory. The relations define the properties behind the shock in terms 
of the shock angle (e) for a given freestream condition. The solution is 
obtained by choosing a shock angle, solving the oblique shock relations, and, 
with this knowledge, finding q^ in Eq. (2.33). Eq,(2.34) must also be satisfied 
and is used to gauge the inaccuracy of the shock angle. An iterative solution 
driving the left-hand side of Eq. (2.33) to zero produces the desired result. 

The distinction between ideal, frozen, or equilibrium solutions comes 
into play only in the oblique shock solution and the relationship between q 
and A. 

To check the accuracy of the approximate treatment, a comparison was 
made against the conical flow tables given in Ref. 20, An examination of the 
results presented in Fig, 2“2 reveals excellent agreement over the entire range 
of approach Mach numbers and cone semi-vertex angles. 
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Fig. 2-2 - Variation of Mach Number at the Surface of a Cone with Cone Semivertex Angle 
for Various Upstream Mach Numbers 
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Section 3 
RESULTS 

Heating rate data correlations for cylinders and cones in laminar 
and turbulent flow are presented in this section. Also presented for com- 
parison with the data are results of the theoretical investigations conducted 
to provide a meaningful extrapolation of the data to flight conditions, 

3,1 YAWED AND UNYAWED CYLINDERS 

Yawed cylinder data correlations were made using the following four 
methods, i.e,. 


Nu^ 

versus 



versus 


Nu 

00 

versus 

Re 

00 

St 

00 

versus 

Re 

00 


From the observed trends of the data, there was no advantage seen in corre- 
lating the data in forms other than Nu versus Re ; therefore, this is the 

00 00 

only form presented. Cylinder test data were obtained from numerous 
sources for laminar flow. However, for turbulent flow, very little test data 
were available. 
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3.1.1 Laminar Results 

The resulting data correlation for laminar flow over cylinders at zero 
degree angle of attack is shown in Fig. 3-1, Sources of the data are noted 
on the figure. As the figure shows, the data correlate well with the Blasuis 
flat plate results, namely 


Nu 


0.332 


V 


Re 


Pr 


0.4 


Figure 3-2 shows the Stanton number versus Reynolds number correlation 
for the same data. Fairly good agreement is also noted with the Blasius 
flat plate results. Figures 3-3 through 3-6 show the freestream Nusselt 
number versus freestream Reynolds number correlations for angles of attack 
of 10, 20, 30, 40 and 50 degrees. Again the sources of the selected data are 
noted on the figures. Also shown are two yawed infinite cylinder theories. 
Beckwith’ s laminar cylinder theory. Ref. 16 is shown as the dashed line, 
whereas the data correlation theory of Ref. 15-Kemp, Rose and Detra, is 
shown as the solid line. The theory of Kemp, Rose and Detra is shown to 
correlate best at all angles of attack except 50 degrees, where the difference 
between the tv^^o theories is small. Referring to the figures, Mach number 
effects are evivdent. 

3.1.2 Turbulent Results 

As mentioned earlier, the available test data for turbulent flow over 
yawed cylinders was very limited. Figure 3-7 shows the Nusselt number 
versus Reynolds number data correlation for angles of attack of 30 and 50 
degrees. As is evident from the figure, Beckwith’s turbulent infinite cylinder 
theory correlates well with the test data. 
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3.2 YAWED AND UNYAWED CONES 

Cone test data correlations were made using local and freestream flow 
properties. The local flow properties were calculated using the equivalent 
cone method. Cone surface pressures predicted using this method for the 
stagnation line were compared with available test data, and in all cases good 
agreement was evident. L»ess data scatter was evident using the local flow 
properties, therefore only these results are presented. 

3,2.1 haminar Results 

Figure 3-8 shows the local Nusselt number versus local Reynolds 
number data correlation for several cone half angles and Mach numbers, 
all at zero- degree angle of attack. Note that most of the data are grouped 
together in a line and correlated well when plotted in this form. The wall 
temperature-tO“total-temperature ratios for the data ranged from 0,2 to 0,5, 
The Mach number effects (for the tunnel conditions used) were small. The 
Mach number varies from 2.49 to 10.1. The cone half angle varied from 5 
to 25 degrees. However, from a plot su Ji as this, where more than one 
parameter is changed between sets of data, it is difficult to isolate the indi- 
vidual effects. In Section 4 of this report it is shown that Mach number and 
cone half angle effects are present and should be considered. 

Also shown in Fig. 3-8 are the theoretical results obtained using Eqs. 
(2,1) and (2,2), There is very little difference between the two theories 
although Eq. (2.1) tended to predict the data slightly better than Eq. (2.2). 
Figure 3-9 shows the local Stanton number versus local Reynolds number 
data correlation. Figures 3-10 through 3-20 shows the data correlations for 
cones in laminar flow at angles of attack from 5 to 45 degrees. Figures 
3-10, 3-12, 3-14, 3-16, 3-18 and 3-20 show that the trend of the data with 
angle of attack is predicted quite well using Streamline Divergence theory. 
However, to obtain the theoretical variations of heating rates with angle of 
attack, Eq, (2,4) was not used per se. Using Eq. (2.4) per se, the theo- 
retical results showed a huge increase in heating rate as the angle of attack 
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increased from zero degrees. This resulted in theoretical results approx- 
imately 70% higher than the data at all angles of attack. However, the theo- 
retical results showed that the proper variation with angle of attack could be 
predicted by Eq. (2,4) if this initial discrepancy could be corrected. For 
this reason the following equation for was used and showed good cor re - 

i-f 

lations with the data, i.e., 


X 


^ a a = 0 deg 


l + 2(l + NEj^)^ 




( 2 . 35 ) 


where 


X. 


"a = 0 deg 



The data of Bushnell (Ref. 3) tend to fall below the theoretical results at all 
angles of attack. This is attributed to the data being taken on a blunt cone, 
and the method used (in Ref. 3) to obtain the local flow properties for use in 
reducing the heat transfer data was different than that used in this report. 

An estimation made to determine the magnitude of the flow property effects 
showed that the data should be about 15% lower than the theoretical sharp 
cone results. 

3.2.2 Turbulent Results 

Figures 3-21 through 3-26 show the turbulent cone data correlations. 
Turbulent cone data available in the literature were more scarce than laminar 
cone data. Shown in Fig, 3-21 along with the data correlations are the theo- 
retical results obtained using Eqs. (2.10) and (2,13), The theoretical results 
obtained using Eq. (2.11) fell between the results obtained using Eqs. (2.10) 
and (2.13) and in most cases very small differences were observed between 
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the three theories for the flight conditions investigated. Figure 3-21 shows 
that the results obtained using Eq. (2.10) tend to correlate best with the test 
data. It can be seen that the Mach number effects are well predicted by the 
theories and they are more pronounced for turbulent flow than for laminar 
flow. The freeflight data of Ref, 7 was predicted well by Eq. (2.10) through- 
out the portions of the flight envelope where the flow was completely turbu- 
lent. From Figs. 3-23 and 3-25 it is evident that the Streamline Divergence 
theory, Eq, (2.17) predicts the data trend quite well at angle of attack. No 
problems were encountered using Eq. (2.17) per se as were encountered in 
the laminar case. 
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a = 0 deg 

O LMSC/805538, Hemispherical Nose Configuration, = 10 (Ref.ll) 

0 LMSC/805538, Flat-Nose Configuration, M =10 (Ref. 11 ) 

/\ LMSC/A838943, Agena Configuration No. 2, M^= 6, 8, 10 (Ref. 1) 

^ LMSC/A838943, Agena Configuration No, 1, M =6,8, 10 (Ref. 1) 

— Flat Plate Theory 


Freestream Nusselt Number vs Freestream Reynolds Number for Laminar 
Flow over Cylinder at a = 20 Degrees 
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NOTES: 
a = 10 deg 

© L»MSC/805538, Hemispherical Nose Configuration, = 1C 

LiMSC/805538, Flat Nose Configuration, = 10 
A LMSC/A838943, Agena Configuration No. 1, M =5, 6, 8, 10 


Fig. 3-3 - 


Free stream Nusselt Number vs Free stream Reynolds Number 
for Laminar Flow over Cylinders at a = 10 Degrees 


LMSC/HREC D162315 













NOTES 


L.MSC/805538, Hemispherical Nose Configuration, M =10 
•LMSC/805538, Flat Nose Configuration, M =10 
LiMSC/A 838943, Agena Configuration No. 1, =3, 5, 6, 8, 10 

LMSC/ A838943, Agena Configuration No. 2, M =3, 5, 6, 8, 10 
Theory, Beckwith, NASA TR-R-40 
Theory, Rose, Detra, Kemp 


Fig. 3-4 - Free stream Nusselt Numbe 
for Laminar Flow over Cvli 
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Freestream Nusselt Number vs Freestream Reynolds Numbei: 
for Laminar Flow over Cylinders at a = 30 Degrees 
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Freestream Nusselt Number vs Freestream Reynolds Number 
for Laminar Flow Over Cylinders at a = 40 and 50 Degrees 
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Local Nusselt Number vs Local Reynolds Number for Laminar Flow 
over Cones at a = 0 Degrees 
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Local Nusselt Number vs Local Reynolds Number for Laminar* Flow 
over Cones at a = 5 Degrees 
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Fig. 3-15 - Local Stanton Number vs Local Reynolds Number for Laminar 
Flow over Cones at a = 15 Degrees 
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Local Nusselt Number vs Local Reynolds Number for Laminar 
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NOTES 


20 deg 
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Local Stanton Number vs Local Reynolds Number for Lamina 
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NASA TN D-4792. 0 = 
— Theory H*/Blasius, 0 


Local Nusselt Number vs Local Reynolds Number for Lamihar 
Flow over Cones at a = 30 Degrees 
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Theory H*/Blasius 


Local Stanton Number vs Local Reynolds Number for Laminar 
Flow over Cones at a = 30 Degrees 


BBSSBa 

B 

5 

sss 

B 

B 


IB 


LMSC/HREC D162315 











Local Nusselt Niimber vs Local Reynolds Number for Laminar 
Flow over Cones at a =45 Degrees 
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Local Nusselt Number ve Local Reynolds Number for Turbulent 
Flow over Cones at a = 0 Degrees 
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NASATND-179, 6 =10 deg, M =4.95 i 

NASA MEMO 10-8-58L, 0 = 10 deg, M^ = 1.98 ■ 

NACA MEMO RM L57D04, 0 = 25 deg, M^ = 1.67,3.97 (Freeflight Datay' 
NASA TN D-824, 0 = 7 l/2 deg, M = 2.8, 3.09 (Freeflight Data) ■ 


Local Stanton Number vs Loca^,Reyno^ds Number for Turbulent Flow 
over Cones at a = 0 Degrees 
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Local Nusselt Number vs Local Reynolds Number for Turbulent 
Flow over Cones at a = 10 Degrees 
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Local Stanton Nvunber vs Local Reynolds Number for Turbulent 
Flow over Cones at a = 10 Degrees 
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20 deg 

NASA TN D-179, 0 = ; 
— Theory H*/Blasius, 0 


Local Nusselt Number vs Local Reynolds Number for Turbulent 
Flow over Cones at a = 20 Degrees 
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Section 4 

CONVECTIVE HEAT TRANSFER DESIGN CURVES FOR CONES 
AND CYLINDERS AT ANGLE OF ATTACK 

The purpose of this section is to present the results of the work of 
Sections 2 and 3 in a form that can be used to determine heating rates to cone 
and cylinder shapes in f reef light at angle of attack. From the design curves 
presented, it is possible to obtain windward streamline and off -windward 
streamline heating for both laminar and turbulent flow. The following sec- 
tions discuss and explain these design curves. Then, two example problems 
are given, one for a cylinder and one for a cone. Following this, a discussion 
of the range of applicability of these "isolated*' cone and cylinder curves to 
composite cone -cylinder-flare shapes is given. 

4.1 HEAT TRANSFER DESIGN CURVES 

For the convenience of the user these curves are presented in the 
following order: 


Type Curves 


Figure Numbers 


I. Cylinders 

A. a = 25 deg, laminar and turbulent, M = 2 to 12 

00 

B. Angle of attack correction factors for a / 25 deg 
(except a = 0 deg) 

C. a = 0 deg, laminar and turbulent 

D. Circumferential distributions, laminar and 
turbulent 


4-1, 4-2 

4-3, 4-4 
4-5 

4-6, 4-7 


II. Cones 

A. a = 0 deg, laminar, M = 2 to 12 

00 

B. a = 0 deg, turbulent, M = 2 to 12 

00 


4-8-4-11 

4-12-4-16 
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Type Curves Figure Numbers 

II. Cones (Continued) 

C. Angle of attack correction factors for a / 0 deg, 

M = 2 to 12, laminar and turbulent 
00 

D. Circumferential distribution, laminar and 
turbulent 

III. Flowfield Parameters 

A. Cylinders, = 1.5 to 26, a = 5 to 50 deg (altitude 
= 100,000 ft) 

B. Cylinders, altitude correction factors 

C. vCones, = 1.2 to 26, a + 0 = 5 to 50 deg (altitude 
= 100,000 ft) 

D. Cones, altitude correction factors 

E. Viscosity and conductivity of air as a function of 
temperature and pressure 

F. Adiabatic wall temperature 

IV. Range of Applicability 

The details of each of these curves are now discussed beginning with 
the cylinders. Freestream Nusselt number, at a = 25 deg, as a function of free- 
stream Reynolds number is presented for laminar and turbulent flow on Figs. 
4-1 and 4-2. It is noted that there is a Mach number effect of these values, 
therefore, separate curves are shown for = 2 to 12. These curves are 
based on a wall temperature of 560^R, This wall temperature affects the 
Nud ^ value only through the heat transfer coefficient in the Nu definition 

(h^D/k^). Note that the k value used here is the freestream value making 

Nu a freestream Nusselt number. An angle of attack value of 25 deg was 

chosen as a ** reference" condition for presenting these curves. For angle 

of attack values other than 25 deg, the curves of Figs. 4-3 and. 4-4 are 

provided — again for laminar and turbulent flow. The curves give the ratio 

of Nup ^ at a given angle of attack to NUj^ ^ at an angle of attack of 25 deg. 

Since there is a Mach number effect, curves are presented for M = 2 to 12, 

00 

4-2 
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4-17, 4-18 
4-19, 4-20 


4-21 - 4-23 
4-24-4-25 

4-26-4-28 
4-29, 4-30 

4-31, 4-32 
4-33 
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These curves are also presented for a wall temperature of 560°R, Figures 4-3 
and 4-4 are for angles of attack down to 5 deg because the Nu based on diameter 
does not correlate at a values near 0 deg. The correlating parameter near 
0 deg is Nu based on x, running length down the cylinder, versus Reynolds 
number also based on x. Therefore, Fig. 4-5 is given for use at a = 0 deg. 

Both laminar and turbulent curves are provided. 

Circumferential heating distributions on cylinders at angle of attack 
values of 10 to 50 deg are given for laminar and turbulent flow in Figs. 4-6 
and 4-7. Specifically, heat transfer coefficient at an angle measured away 
from the windward streamline to the local point, ratioed to the windward 
streamline value is shown as a function of These figures were derived 
from the data of Ref. 1 for laminar flow, and from Ref. 8 for turbulent flow. 

The cone heating curves are shown beginning with Fig. 4-8. The heating 
parameters for cones^ are correlated and presented in a way similar to that 
used for cylinders. However, there is one additional parameter to be con- 
sidered, namely the cone angle. Also local rather than freestream properties 
are used, and the reference*' value of Nu is at a = 0 deg rather than 25 deg. 
Local Nusseit number at a = 0 deg based on x , the distance from the cone apex 
to the local point at which heating rates are needed, versus local Reynolds num- 
ber are shown on Figs. 4-8 to 4-16. Figures 4-8 through 4-11 are for laminar 
flow while 4-12 through 4-16 are for turbulent flow. Note that a separate figure 
is given for each Mach number and that varying cone angles are shown. 

For cones at angles of attack other than 0 deg, the correction factors of 
Figs, 4-17 and 4-18 are given for laminar and turbulent flow, respectively. 
These curves present the ratio of local Nusseit number at a given angle of 
attack to the local Nusseit number at an angle of attack of 0 deg as a fu,nction 
of angle of attack. These values are given for cone angles of 5, 15 and 25 deg 
and for Mach numbers of 2 to 12. 
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These correction factors were generated for a specific altitude of 
100,000 ft, however, the altitude correction factor is small for altitudes from 
sea level to approximately 400,000 ft. These curves appear to be somewhat 
random which is due to two factors. First is the fact that the shock angle on 
the cones varies with three parameters, cone angle, Mach number, and angle 
of attack. The second factor is the erratic behavior of thermal conductivity 
(Fig. 4-32) with temperature and pressure in the regions covered. 

Circumferential heating distributions for cones are presented in Figs, 
4-19 and 4-20 for laminar and turbulent flow. As in the case for cylinders, 
the heat transfer coefficient ratios are presented as a function of <f> , the angle 
measured away from the windward streamline. These curves are given as 
reasonable approximations for use for cone angles up to approximately 25 deg 
rather than for any specific cone angle value. A scarcity of cone data did not 
allow a comprehensive determination of the effect of cone angle during the 
present study. The results shown in Figs. 4-19 and 4-20 are fairings and ex- 
trapolations of the data in Refs. 3 and 9 for laminar and turbulent flow, 
respectively. 

4.2 LOCAL FLOWFIELD CURVES FOR CONES AND CYLINDERS 

In order to obtain Nusselt numbers, hence heating rates, from the pre- 
ceding design curves it is necessary to know the local Reynolds number. 

To assist the user in obtaining these, the curves of this section are presented. 
First, for cylinders, local static to freestream static temperature, density, 
and velocity ratios are given in Figs. 4-21, 4-22 and 4-23. These are given 
as a function of Mach number and angle of attack. Note on Fig. 4-23 that 
the velocity ratio variation with Mach number is negligible. These figures 
are for an altitude of 100,000 ft. For different altitudes Figs. 4-24 and 4-25 
present correction factors as indicated on each respective curve. 

For flow fields on cones. Figs. 4-26, 4-27 and 4-28 give temperature, 
dem|ity, and velocity ratios as a function of Mach number. Here the sum of 
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angle of attack and cone half-angle has been used rather than angle of attack 
as in the cylinder curves. Again, these curves are for an altitude of 100,000 
ft, and altitude correction factors are given in Figs. 4-29 and 4^30. 


Figures 4-31 and 4-32 present viscosity and thermal conductivity for air 
as a function of temperature and pressure. Viscosity, of course, is needed 
in determining Reynolds number, and conductivity is needed for determining 
heat-transfer coefficient from the Nusselt number. Figure 4-33 presents the 
ratio of adiabatic wall temperature to total temperature, as a function of free- 
stream velocity. The total temperature used here, is defined to be the ideal 
gas value as follows: 


V 


00 


^ 2g J G 


P 


This is the value used to obtain to be used in conjunction with h^ 
obtained from the Nu to get heating rates. This is illustrated in the following 
example problems, 

4.3 APPLICATION OF DESIGN CURVES 
4.3.1 Example Problem No. 1 

To illustrate the use of these curves the following example problem is 
given. Assume a cylinder with a 2-ft diameter flying at the following conditions 


Altitude = 150,000. ft 
Velocity = 13,000. ft/sec 
Angle of attack = 40. deg 
W all tempo ratu re = 5 60° R. 
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Heating rates are desired on the windward streamline and 45 deg away from 
the windward streamline. First, the frees tream Mach and Reynolds number 
must be determined. 

From the 1962 Standard Atmosphere table, the following are determined 


M 

00 

= 12.1 

Per. 

00 

= 1.11 X lO" 

T 

00 

= 479. 0°R 

p 

00 

= 1.34 X 10" 

L jtx 

= l.lx 10’® 


m' 


3.6 X 10 ^ Btu/ft-sec-^R. 


m' 


00 


From these values 


P V D , 

= 2.62 x 10' 


At this Reynolds number, the flow should be laminar, therefore from Fig, 4-1 


Then from Fig. 4 --3, 


Nu~ = 8.6 X 10 . 

D, 00 , a = 25 deg 


Nu 

Nu 


p, ao , a = 40 deg _ j ^ 
D, do, a = 25 deg 
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and 


Nu„ . = 8.6 X 10^ X 1.7 = 14.6 X 10^ 

D,oo,a=40deg 


The neat transfer coefficient is now found from 


_ ^00 _ 14.6 X 10 X 3.6 X 10 _ ^ 1 rt-3 Btu 

— = — — ■ ■■ = 2.63 X 1 0 S ' ■ ■ ■ ■ ■ 

c D 2 • ,, 2 o. 


ft -sec- R 


From Fig. 4-33 (at = 0.85 for laminar flow) 

T 


T = 12,700^R 
aw * 


and from 


q = h (T -T ) 
^ c aw w' 


-3 


q = 2.63i X 10'"^ (12,700.- 


Q = 32.0 


-sec 


This is the heating rate on the windward streamline. 


For the point at 45 deg away from the windward streamline. Fig. 4-6, is 
used, where 


h. = 45 deg 


= 0 deg 


= 0. 


Therefore, 




^ = 45 deg 


= 0.66(32.0) = 21.1 


Btu 


ft -sec 
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= 0.99 for 50,000. ft 

= 1.01 for 50,000 ft 


Therefore, for 50,000 ft 


and 


~ = 0.99(8.2) = 8.12 
^1 


= 1.01 (6.2) = 6.46 


T, = 8.12(390.0) = 3160.0°R 


p2 = 6.46(1.16x10“^) 


= 7.28x 10 


, lb 
-2 m 


ft' 


V 2 = 0.785 ( 1 0,000) = 7 , 850 ft/sec 


-5 


From Fig. 4-31, ^2 = 3-9 x 10 


From Fig. 4-32, = 1.5 x 10 


-5 


Btu 


ft-sec- R 


From which 


Re, T = ? « 14.6 x10^ 

x,L 2 
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At this Reynolds number it can be safely assumed that turbulent flow exists. 
From Fig. 4-16, 


Nu 


X, L 


i-' 


: 4,2x10 

0 deg 


M = 10 deg 
00 


For an angle of attack of 20 deg and M = 10 

00 


Nu 

rST 


x,L,(y-20 deg from Fig. 4-18 

X, li, a = 0 deg 


and 


«»=<,I,.« = 20d.6 = * ">*• 


From this, the following results: 


■ ■■ ■.* A , 


2 X 10 


-1 


X 


From Fig. (4-33) 


= SOOO.^R at = 0.9 

* - "t ij, 


and 


q = h (T -T ) = 890. 
^ c ' aw w' 


Btu 


ft -sec 


on the windward streamline at a point 10 ft from the apex. 
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Now for the heating at a point also 10 ft from the apex but at 30 deg 
away from the windward streamline* refer to Fig, 4-20. At ^ = 30 deg* a = 
20 deg 

^ =30 deg _ J 
*^0=0 deg 

from which 



= 30 deg 


1.11(890) = 980, -Sill — 

fr-sec 


4.4 REGIONS OF APPLICABILITY OF ISOLATED BODY DESIGN CURVES 

When using the preceding design curves on composite shapes such as 
a cone -cylinder-flare, it is desirable to know how close to a junction the 
curves apply. Two cases are considered here, downstream of an expansion 
corner, and upstream of a compression corner. Data for these regions are 
limited in that instrumentation locations are usually too widely spaced in the 
corner region to give a distribution. However, the curves of Fig. 4-34 are 
given as an approximate means for determining the distance downstream of J 
an expansion corner at which the cylinder curves become applicable. 

Figure 4-34 presents x/d, the distance downstream- to- cylinder diameter 
ratio, as a function of freestream Mach number for various angles of attack. 
This curve is for laminar flow only. However, indications are that when there 
is turbulent flow ahead of a corner, it goes laminar downstream in the expan- 
sion region. Therefore, this curve should give a reasonable approximation. 

On the cone upstream of the expansion corner^ the isolate d^^c^ 
should apply all the way to the corner. 

In compression corner regions, the available heating data were not con- 
sistent enough to allow the development of a set of curves similar to Fig. 4-34, 
In lieu of this, the following equations from Ref . 21 are offered as a nie^s for 

// 
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handling problems in these regions. Upstream of a compression corner, the 
isolated cylinder curves apply up to the separation point. The problem is 
then to define where separation will occur. This can be estimated by defining 
a critical pressure coefficient at which separation will occur. The relation- 
ships are: 

For laminar flow: 


CP 

critical 


2.03 {M 2 - 1) 


- 0.306 




0.25 


For turbulent flow: 


CP = 

critical " /«_ \0.1 




Separation should not occur if CP is less than ^^^ritical defined 

by 


CP = 


P, -P 
d u 

i p 


where P^ is the static pressure downstream of the corner, and P^ is the 
static pressure upstream of the corner. If separation occurs, the following 
expressions can be used to estimate the location of separation: 


fp 1 ^-056 
570.0 ^ 

.1 i" - 
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For laminar flow; 


d 

u 
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where is the distance downstream of the corner to the reattachment point 
For turbulent flow: 

Id k%. „) ''• 


1.2 


u 




3.1T5 


To estimate the increase in heating at the point of reattachment, the 
following relations are used: 







Notes 


1. Flow is laminar 

2. Wall temperature 


Reference Free stream Nusselt Numbei vs Free stream Reynolds Number 
for Laminar Flow over Cylinders at a =25 Degrees 


r i. 


. j. .+ 


: 

-J 

IT - ;■ 

1 * 

^ t 

- i 

rj . :i 


r T 

} i 

!;■■■ r- 

•• L — • -- 4. ^ 


•i 1 
! 

1 1 , : 1 

yLi 

« • 

. i 
• 

: i 


■' j 


1 •.'••-1 

P-TJ 

■ — i 

1 1 

hr- 

h- 

Ll 

r- 

!“1 

h 

r-i 

M 

j 


! ! 






1 

i 

i- 

i ' ' 

1 "1 

1 1 


1 . • 


. • . i 

r-TT- -■ 


( 






r "i 
[ 1 


I 


1: -j 


■ 



\ 

1 .: 

...*i 

k 



ji;:- 

I- 


[0 

i 




1 f j 

1 • • . i 

i ■ 

• 

j\ 

’i'l 
» • * 


i -] 

i -1 
1 ' 1 

‘ 1 

I7 

I 







!. .. 


1 

1 ' 

1 , 

[. 1 

r I 
1 

' 1^ 

; i 



i ■ 

• : • 

■■Mlj 

1 ■' 

r * 1 

1 •• 

1 . *1 

L ' 

L;l 



1 ' 

1 1 

' I 
I. J 

1 1 


f . ■ 1 

1 1 

J 



L-'. ! 


!-Pi 

1 t * ; i 

L 



I I 

I 

1 i 

L_] 

i: 


y . . 1 












Notes 


1, Flow is turbulent 


2, Wall temperature 


Reference Free stream Nusselt Number vs Free stream Reynolds Number 
for Turbulent Flow over Cylinders at a =25 Degrees 
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NOTE 


1. Flow is laminar 

2. Wall temperature 


Angle of Attack, a (deg) 


Angle of Attack Correction Factor for Laminar 
Flow over Cylinders 
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1, Flow is turbulent 

2. Wall temperature = 560°R 


Angle of Attack, a (deg) 


Angle of Attack Correction Factor for Turbulent Flow over 
Cylinders 








£ 



Wall Temperature = 560 R 


ream Nusselt Number vs Freestream Reynolds Number 


minar and Turbulent Flow over Cylinder at a = 0 Degrees 
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NOTE 


fficient Distribution 
pt Cylinder 


Fig. 4-6 - Circumfe 
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Notes 


Lv is laminar 


Local Nusselt Number vs Local Reynolds Number for Laminar Flow 
over Cones at a = 0 Degrees, M =2 



. 





• ; . I 

i : : 


. • i 




; i . ; 

! 

! ■' i 


•\ 

i - 
: ' • 

• ; ! ; 
! I * • 

-'i 

■ i ; '! 



1 

4 

i 



’ 4 

. . • i 

; r i 


:-;l7 ‘L 



-- 

Iv: 

1 


--4 


' 

. 

•1 

I — 

• 

' ■ 
T' 

::t 

1 

1-- 



L:| 


_ 




.* • 


Li- 



: 1 
• 

' ' 


u 

-• 

-j 

4' 

r 

1 

. t 
1 



4-1 — j— 


-1 -I 

'^*!i 


1 - 

1 • 

i'l ■ 
k-'4— 
• • 

; ! 
1 

L' = ’ 

f- — 

K : * f - 

; : • 

rn 

• j 

rr 

— 



_L 

4 
- 1 











= 0 deg 



4-^; 1. a - 0 deg 

I . 

2. Flow is laminar 
]~1 3. V/all temperature 


Local Nusselt Number vs Local Royrolds Number for Laniinar 
Flow over Cones at a = 0 Degrees, Mqq = 4 
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Notes 


{f 2. Flow is laminar 
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3, Wall temperature 


Local Nusselt Numbe 
over' Cones a't a = 0 E 


^ 8 Local Reynolds Number for Laminar FI 

degrees, M - 6 ‘ ' 


! 

E3 

’ • • • 
-Tt 

!• :=J 

il::. 

mam 

1 i ' ’ i 

1 • ‘ 

U‘ 


:xi 





: ' : ' • 
: 1 


♦ 

• 

! i i 

i 

; • 1 t 


1 

1 1 • • • f : : 

1 ! -:L_4 

LJ 


ri- 


» — “ — 
1 1 ; 

j : • 

} : ... 
J. . ' . 

r i i 


L- L-L ‘ 

- « - 

. ' 

I ^ 

1 ■ !’■:! 

1 . : 
j ■ ■ 

• 1 
1 . . 1 . . 

j ■ • ■ j 1 1 ■ ! 

[■T-l: 

1 . . i 

i' - -I 

1 


LMSC/HREC D162315 








0 deg 



Local Nusselt Number vs Local Reynolds Number for Laminar. Flow 
over Cones at a = 0 Degrees, M = 8 to 12 
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Fig. 4-12 - Local Nusselt Number vs Local Reynolds Number for Turbelent 
Flow over Cones at a = 0 Degrees, M =2 
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Notes 


1. a » 0 deg 

2. Flow is turbulent 

3. Wall temperature 


Local Nusselt Number vs Local Reynolds Number for Turbulent 
Flow over Cones at a = 0 Degrees, M = 4 
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Notes 


2, Flow is turbulent 

3, Wall temperature 


Local Nusselt Number vs Local Reynolds Number for Turbulent Flow 
over Cones at a = 0 Degrees, M =6 
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) Notes 


1. a a 0 deg 

2, Flow is turbulent 


3, Wall temperature 


Local Nusselt Number vs Local Reynolds Number for Turbulent Flow 
over Cones at a = 0 Degrees, M =8 








Turbulent Flow 

Wall Temperature = 560°R 


Local Niisselt Number vs Local Reynolds Number for Turbulent Flow 
over Cones at a = 0 Degrees, M =10 and 12 
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3.0 



a (deg) 



20 30 

a (deg) 


6=5° 

NOTES: 

1. Wall Temperature = 560' 

2. Altitude = 0 to 400,000 ft 

15° 



10 20 30 

a (deg) 


Fig. 4-18 - Angle of Attack Correction Factor for Turbulent Flow 
over Cones forM =2,4,6,8,10,12 
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NOTE 


Fig. 4-19 - Circumferential Heat Transfer Coefficient Distribution for Laminar Flow over Cone 
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Turbulent flow 




Fig. 4-20 - Circumferential Heat Transfer Coefficient Distribution for Turbulent Flow 
over Cones 
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Local to Freestream Density Ratio, p^/p 
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NOTES: 

1. Cylinder Flow 

2. M = 2.0 to 26.0 


The Ratio of Local Velocity on Cylinder to Freestream 
Velocity vs a 
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NOTES 


at Altitude = 100,000 fret 
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Fitj. 4-24 - Altitude and Mach Number Correction Factors for Static Temperature <tatios on Cylinders 
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Fig. 4-26 - The Ratio of Local Static Temperature on Cones to Freestream Static Temperature 


VS M 


LMSC/hREC D162315 





I 

\ 




LMSC/HREC D162315 




Sum 

■■■■■ 


NOTE 


For flow over cones 
Altitude = 0 to 400, 000 ft 


The Ratio of Local Veloc 

















<1 


ess 


'='H»aniii» 


|gH!lili9i!ili!!»iilL\l 

inyiiiyilSiyiilliiliiigi' 

liriliLiiiiiiilPiOSiaK' 


I^StiiSiSISM 


NOTES 


at Altitude 


at Altitude 


For Flow on Cones 


Altitude (ft) 


200,000 


Altitude (ft) 


Altitude (it) 


Altitude and Mach Number Correction Factors for Static Temperature Ratios on Cones 









CF 





Viscosity (Ib^/ft-sec) 



LMSC/HREC D162316 














Adiabatic Wall Temperature (1000°R) 




8 12 16 20 24 
V^, Freestream Velocity - 1000 Ft/sec 

Fig. 4-33 - Adiabatic Wall Temperature vs Freestream Velocity 
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Section 5 
CONCLUSIONS 

The cylinder test data correlated well when plotted in the form of free- 
stream Nusselt number versus freestream Reynolds number. The effects of 
Mach number were considerable and were predicted well by the theoi^. Lam- 
inar convective heating rates to cylinders at zero degree angle of attack were 
well-predicted by flat plate theory. For cylinders at angle of attack, yawed 
infinite cylinder theory applied. The laminar stagnation line heating rates 
were well predicted by the method of Kemp, Rose and Detra (Ref. 15). The 
turbulent stagnation line heating rates were well predicted by the method of 
Beckwith and Gallagher, Ref. 8. 

The cone test data correlated best when plotted in the form of local 
Nusselt number versus local Reynolds number. The effects of Mach number 
were more pronounced for turbulent flow than for laminar flow. The heating 
rate results were well predicted by Streamline Divergence theory. 

Much more useful test data are needed especially for turbulent flow 
over yawed cylinders and cones. Also, more test data are needed on circum- 
ferential heating rate distributions for a wide range of Mach numbers and 
cone half -angles and for turbulent flow over cylinders. The use of a consis- 
tent method for obtaining local flow properties on cones helped considerably 
in obtaining good data correlations. Extrapolation of ground test data to 
flight conditions is a crucial step. Generally ground test data are obtained at 
low stagnation temperatures, such that real gas effects are not important. It 
is felt that the use of theories for extrapolation of data to flight conditions is 
meaningful.^^^ 
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